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We consider the parity invariant QED2+1 where the matter is represented as a mixture of fermions 
with opposite spins. It is argued that the perturbative ground state of the system is unstable with 
respect to the formation of magnetized ground state. Carrying out the finite temperature analysis 
we show that the magnetic instability disappears in the high temperature regime. 



o 



(3 . I. INTRODUCTION 

o . 

Interest to the (2-|-l)-dimensional field theory is heated up by its successful applications to the problems of planar 
condensed matter systems like the quantum Hall fluid and unconventional superconductors. One of the relevant issues 
concerns the phenomenon of spontaneous magnetization in the systems of planar fermion matter. 

Different authors have given convincing arguments that the perturbative ground state (the one with F^j^u = 0) 
of a 2D dense fermion matter exhibits the magnetic instability. In other terms there does exist the true ground state 
with the nonzero value of a proper magnetic field. 
T-H , The peculiar property of planar physics is that spin is not quantized, and spin-up and spin-down states belong to 
' the different representations of the rotation group 5*0(2). Remark, that the existence of magnetized ground state 
\l ' was demonstrated in the systems of fermions with one and the same spin orientation. Hence one can propose the 
simple mechanism explaining this phenomenon: matter consisting entirely of either up or down spin particles possesses 
I nonvanishing spin density leading to a proper magnetization of the system. 

. In the present paper we study the double-spin model (also referred to as duplicated) where the matter represents 
the mixture of opposite spin fermions. Such a consideration can be motivated by the observation'' that the magnetic 
: properties of duplicated and single-spin (up or down) systems can be quite different. 

' In the rest part of this section we remind some principal aspects of spontaneous magnetization in the single-spin 
-4^ ] system. In Section 2 we perform the perturbative study of duplicated system and find out magnetic instability of 
ft) ^ perturbative ground state. The mechanism of instability in the double-spin model turns out to be different from that 

I in the single-spin systems. Some necessary calculations are placed in the Appendix. 
<^ . The basic object of our account is the finite temperature effective action W[A\ which arises after integrating out 
ILJ ' the fermion fields. In the Matsubara formalism it is given by 



X 
J3 



exp{-W^[yl]} = j Dil^*DiP expj-y idrdrj , 



where L is the Euclidean Lagrangian, < r < T"^, and T is the temperature. 

Consider the nonrelativistic fermion matter interacting with U(\) gauge field. The corresponding Euclidean La- 
grangian is given by 



+ — \{dk + ieAk)4^\^ - Bri^ - i^ri^. 

Im Im 

where p is the neutralizing charge density, a describes the spin degree of freedom, magnetic field is defined as 
B = d2Ai — diA2, and p is the chemical potential. 
The linear part of W[A] appears as 

Wi = e [{iWi^) - p]Ar- ^ {r^)B}dTdr, 
J 2m 

where the symbol (• • •) denotes the thermal average corresponding to F^i, = 0. 
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From the last expression one can extract the principal conclusion concerning the magnetic properties of the single- 
spin system. Due to the term, linear in the magnetic field, the extremum of W[A] will be realized for some B ^Q. The 
common conclusion achieved in the finite temperature analysis of single-spin models, implies that the spontaneous 
magnetization persists at any temperatures. The reason is quite simple: the finite spin density {a /2){il>*^) creates 
the proper magnetization of the system 

ere , ,^ ,. 

which survives at any temperatures, provided the particle density is finite. 

Concerning the ylT--term, it leads to the electric field corresponding to a local charge distribution e{tjj*ip) — ep. 

II. DOUBLE-SPIN SYSTEM 

In this section we consider the duplicated model where the matter represents the mixture of opposite spin (icr) 
fcrmions. 

The system is assumed to be in contact with a particle reservoir which keeps the total particle density fixed and 
guarantees the chemical equilibrium between spin-up and spin-down subsystems. This can be realized by equalizing 
the corresponding chemical potentials. Then (ip'^ipi) = (V'lV'x): ^^id the linear part of W appears as 



where the term linear in the magnetic field is absent since the duplicated model is parity invariant. Q Therefore, the 
configuration with B = Q represents the extremum of the corresponding effective action. The prior question which 
arises in this connection is the one concerning the stability of this extremum. This issue can be studied in terms 
of the two-point functions calculated within the Gaussian approximation. Considering these functions as integral 
operators we can introduce the corresponding eigenvalues which characterize the spectrum of gauge field fluctuations. 
Stability of the perturbative ground state requires the non-negative definiteness of this eigenvalue spectrum. If the 
spectrum contains at least one negative eigenvalue, then one is faced with an instability leading to the formation of 
the nontrivial ground state. In what follows we adopt this criterion and search for the negative eigenvalues of the 
two-point functions determining the Gaussian part. The later is given by 

f f c'^TTl f 

W2 ^ - F^^Ff^^dx H — / A^{xi)G^^{xi - X2)At,{x2)dxidx2, 

where dx — drdr, and G^^^^x) are the current correlation fmictions. In the Fourier representation they look as 

__ fdk^ f{pE+) - J{I3EJ) 
" J 2^2 (kk')2 + m^e ^ 

yjri ^ IT ^tt ? 

^2 1 2 

Gij — Sij T-^ GrT + — 2 ~ kikj) G , 



^Magnetic field defined in 2+1 dimensions is a pseudoscalar. Therefore, the spin interaction leads to a parity violation in 
single-spin systems. Provided spin-up and spin-down fermions are interchanged under the parity transformation the duplicated 
model is parity invariant. 
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where ^ = 2-KTn with n = 0, ±1, ±2 . . . , and (3 = T ^. The quantities E± and f{z) are given by 



■''(^) - 1 + e--M/T • 

The induced Chern-Simons terms corresponding to the spin-up and spin-down fermions are mutually canceled out 

in accordance with the parity invariance. 

From the definition of E± it follows that E± > for sgn(kk') = ±1, and since f{z) is a decreasing function, we 
get 

[f{l3E+) - f{pE_)]{W) <Q. (2) 

This implies that Grr > 0, and the corresponding contributions to W2 are positive. The rest part of W2 is related 
with pure magnetic configurations and determines the two-point function 

S^W2 

dB(xi)dB(x2) 

which, up to a constant factor, is the susceptibility of the pcrturbativc ground state. 

Because of the translational invariance, it is convenient to analyze the eigenvalues of A(a;i,a;2) in the Fourier 
representation. These eigenvalues can be written in the following form 

A(^,k,r) = 1 + ^G(^,k,T), (3) 

where G(^,k, T) is given by the Eq. (1). 

From these expressions one can detect the possibility of the existence of negative eigenvalues leading to the magnetic 
instability. In fact, due to the relation (2), the cr-dependent part of G is negative, and for sufiiciently large values of 
(7^ the quantity G(^, k, T) will become negative at least for some ^ and k. 

In order, to expose the underlying physics, consider the dense matter of spinless fermions. In the absence of 
gauge fields fermions are organized in a Fermi sphere. Turning on the homogeneous magnetic field, fermions will be 
rearranged into the Landau levels, and the orbital diamagnetism will lead to the increase of the energy of the system. 
Let us now attach the spin to these fermions. Inclusion of the spin with positive (negative) aeB decreases (increases) 
the energy. Due to the chemical equilibrium between the opposite spin subsystems, the partial density of fermions 
with creB > exceeds the one corresponding to aeB < 0. Therefore, the spin effects lead to the spin paramagnetism 
promoting the decrease of the energy. If is large enough, then the paramagnetism dominates over the diamagnetism 
and causes the overall decrease of the energy of the fermion system. 

Further, for sufiiciently large values of e^/m, the negative values of G will dominate in Eq. (3), and for some ^ and 
k we shall get A < 0. 

In the light of these arguments we can formulate two main points required for the existence of negative modes: 

(i) the magnitude of a must be large enough to generate the instability of the fermion matter. 

(ii) the value of the fraction e^/m must be large enough to guarantee that the matter instability will not be 
compensated by the gauge field contributions. 

In the following two subsections we present more detailed analysis of the negative modes for static (^ = 0) and 
nonstatic (^ ^ 0) cases, respectively. 

A. Static case 

Let us first explore the point (i) concerning the negative values of G. As it is shown in the Appendix, the quantity 
G(0, k, T) can be written in the following form 

'^"■''■^'°r?? i+e<»"-.c./r . w 
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where C = k'^/SmT and k = |k|. 

Assuming that the total particle density is fixed, we express the chemical potential in terms of other parameters. 
The defining equation is + (V'JV'j,) = f^e, with 



mT 



and consequently 



T In (e^/® - 1 



where Q = mT/irne. 

Consider first the zero temperature limit. In that case the Fermi distribution function becomes steplike, and we get 



G(0,fc,0) = 



1 - 3^2 



1 



1 



1 - 3cr 



A-2 

A;2 



0{k - kp), 



127r Utt 
where the characteristic momentum kp is defined by 

k% = 8m lim u, = SnUe ■ 

In Fig. 1(a) we present (7(0, k, 0) for the different values of a^. As one can be convinced, the negative eigenvalues 
appear only for > 1/3. 



G 



= 




a2 = 1/6 , 




(72 = 1/3 j 




(72 = 1/2 1 




(72 = 1 
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1 k/kp 1 

FIG. 1. (a) G(0,fc,0) versus k; (b) G(0,A;,T) versus k, cr^ = 1. 

Consider now the point (ii) and trace out the condition guaranteeing that the negative values of G(0, k, 0) will lead 
to A < 0. For > 1/3 the minimal value of G(0, fc,0) is given by 



G(0,0 < k < kp,0) 



3a' -1 



and the required condition appears to be 



m 3(7^-1 
? ^ 127r 



(5) 



In this case one obtains that A(0, k, 0) < for < A; < fcc, where kc is some critical value which is the solution to 
A(0, fc,0) = 0. Remark, that kc > kp. 
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The finite temperature beliaviour of G(0, k, T) is depicted in Fig. 1(b), where the case of cr^ = 1 is represented for 
different values of the parameter 0. We see that G(0, fc,T) tends to zero when T increases. Therefore, the negative 
modes, observed at low temperatures, disappear in the high temperature regime. 

Remark, that the relation (5) can be realized only for > 1/3 and therefore, embraces the points (i) and (ii) 
simultaneously. Thus, the relation (5) appears to be a sufficient condition for the existence of A(0, fc, T) < 0. Moreover, 
it is also a necessary one, since in the opposite case A(0, k, T) is positive for all k and T. In order to check up this 
assertion, one can carry out some estimates. Note, that in (4) we have < s < 1. In this interval e^^"'^"'' ^ > e^^^~^^ 
and e^(^~* ^ < e^^^~* \ Using these inequalities in (4), one gets 



A(0,fc,T) > 1 + — - In- 



Further, since C > 0, we have 



and consequently. 



m Un C l+c/'^-f 



. 1 , 1 +e^/^ 
< - In- — ,^ , < 1 



C l+e'^/'^-C 



A(0,fc,T) > 1+ - 

m 12-K 

This is a general relation valid for all values of k and T, and leading to A(0, fc, T) > when the condition (5) is not 
held. 



B. Nonstatic case 



Integrating out the polar angle in (1), the corresponding expression can be presented in the following form 



G(C,fc,T) 



An 



(^-i)/(C^'), 



(6) 



1 4to2^2 _^ ;,4(i _ ^2) I 

2 4to2^2 + ^4(1 ^_ ^)2 + Ye' 



In Fig. 2 we depict the behaviour of G(^, k, T) for a"^ = 0. In that case paramagnetism is absent and the eigenvalue 
spectrum is positively defined. 
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FIG. 2. (a) G(i,k,T) versus k,(, = 0.1(7rne/m), a'^ = 0; (b) G{^,k,T) versus ^, k = O.Skp, a'^ = 0. 



The case tr^ = 1 is depicted in Fig. 3, where we observe that G becomes negative for some values of ^ and k. 
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FIG. 3. (a) G{^,k,T) versus k, ^ = 0.1(7rne/m), a'^ = 1; (b) G{^,k,T) versus ^, A; = 0.8A;f, ct^ = 1. 

Prom the expression (6) one can extract the asymptotic properties of G{^, k, T) for fc — > oo and ^ — > oo. In the 
large-fc hmit we obtain 



2nd fc2 



for 



G ^ 8 



/ ln(l + e--+''/^)rf^ for 
Jo 



7r52 k^ 



where ^ = 1 + (2m^/fc^)^. These expressions are vaUd irrespectively whether the fraction ^/k"^ vanishes, diverges or 
stays finite when fc ^ oo. 

Consider the limit ^ ^ oo. The case when k tends to infinity together with ^ has been already discussed above. 
Therefore, we assume that k is finite and get 



G 



1 - (7^ TP 



ln(l + e^/^) + - 4t / + e-^+^Z^) 

Jo 



dz . 



The asymptotic relations imply that the negative values of A can be located in a finite region of the (^, fc)-plane. 
Figures 2 and 3 demonstrate that the magnitude of the negative values of G tend to zero as T increases. Therefore, 
in the high temperature regime we get A > 0. 
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IV. APPENDIX 



Here we comment on the calculation of G(0, k, T). The common details are shown for its a-dependent part: 

dk' f{pE+) - f{0E_) 



C(0,k,T) 



f dk' 
1 2^ 



kk' 
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Due to the singularity, the E± parts are separately divergent and the integral cannot be decoupled in the corresponding 
way. 

Expanding f{(3E±) in powers of kk', we use the polar variables and integrating over k' get 



o„(0,t,r)=-i:±^/.».(C) 



n=0 



where C = k'^/SmT. Expand /("^(C) and arrange the powers of ^. Besides, we use /(J)(0) = (-l)J>(j)(/z/T) with 
(l){z) = (1 + e"^)-^ and get 

1 ^ / 1 \n/-n 

G„(O,t.T)=^|:jl^0WWr). (A^l) 

Analyze the structure of (j)^"-\z). From (^'{z) = </>(! — (p) it follows that can be presented as an (n + l)'th 
order polynomial with respect to (j) 

(z) = P„+i[0(z)]. (A.2) 

These polynomials satisfy the recurrency relation P„+i = 0(1 — (j)){dPn/d<j)) with Pi[0] = By the direct calculation 
it can be checked up that the solution to this recurrency chain appears as 

oo k 

P„+i[</-] =5^^(-l)'(Z + irc[/+^ (A.3) 

A;=0 1=0 

Substituting Eqs. (A.2) and (A.3) into Eq. (A.l), we use 



ra=0 

and summing up over I, get 

[l-e-C(i-^^)]fc 
[1 + 6-^/^]'^+^ 



Jo 



The infinite sum over k converges, yielding 

ds 



GrriO, k,T)= I f 



Performing the similar manipulations for the a- independent part of G and combining with (A. 4), we get (4). 
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